aves<-read.table("aves_cerrado.csv", header=T, sep=";")
head(aves)
aves$urubu[is.na(aves$urubu)] <- 0
aves$carcara[is.na(aves$carcara)] <- 0
aves$seriema[is.na(aves$seriema)] <- 0
aves$fisionomia[aves$fisionomia=="ce"] <- "Ce"
aves$fisionomia <- factor(aves$fisionomia, levels=c("CL","CC","Ce"))
str(aves)
head(aves)
#1.Ha relacao entre o numero de avistamentos de seriemas e carcaras em alguma das fisionomias amostradas?
mod.ce<-lm(seriema~carcara, data=aves[aves$fisionomia=="Ce",])
mod.cc<-lm(seriema~carcara, data=aves[aves$fisionomia=="CC",])
mod.cl<-lm(seriema~carcara, data=aves[aves$fisionomia=="CL",])
summary(mod.ce)
summary(mod.cc)
summary(mod.cl)
p.ce<-summary(mod.ce)$coefficients[2,4]
p.cc<-summary(mod.cc)$coefficients[2,4]
p.cl<-summary(mod.cl)$coefficients[2,4]
coef.ce<-coef(mod.ce)
coef.cc<-coef(mod.cc)
coef.cl<-coef(mod.cl)
#Ha diferencas na relacao entre o numero de avistamentos de seriemas e carcaras entre as fisionomias do cerrado?
aves.anc<-lm(seriema~carcara*fisionomia, data=aves)
summary(aves.anc)
#Resposta: Nao ha diferencas na relacao entre o numero de avistamentos de seriemas e carcaras entre as fisionomias do cerrado